Sodium: a charge-transfer insulator at high pressures 
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By means of first-principles methods we analyse the optical response of transparent dense sodium 
as a function of applied pressure. We discover an unusual kind of charge-transfer exciton that pro- 
ceeds from the interstitial distribution of valence electrons repelled away from the ionic cores by the 
Coulomb interaction and the Pauli repulsion. The predicted absorption spectrum shows a strong 
anisotropy with light polarization that just at pressures above the metal-insulator transition mani- 
fests as sodium being optically transparent in one direction but reflective in the other. This result 
provides a key information about the crystal structure of transparent sodium, a new unconventional 
inorganic electride. 

PACS numbers: 



Simple-to- complex phase transition- In textbooks 
sodium is defined as the simplest of the simple metals [l| . 
At ambient pressure it has a high-symmetry crystal struc- 
ture with non directional metallic bondings. Valence 
electrons can be treated as almost independent particles, 
only weakly perturbed by the electron-ion interaction. 
They form a homogeneous and isotropic electron gas, giv- 
ing rise to a nearly spherical Fermi surface. As a conse- 
quence, also the optical properties are well described by 
the Drude model of free carriers. Application of pressure, 
instead, changes the situation quite dramatically while 
revealing an emergent complexity. Sodium under pres- 
sure shows an unexpected and often counter-intuitive be- 
haviour. Instead of becoming more free-electron-like, at 
higher densities it adopts several low symmetry (and even 
incommensurate) crystal structures ^2-7] and has a very 
unusual melting curve 0, S] ■ This surprising behaviour 
had been anticipated by Neaton and Ashcroft 9], who 
also predicted a possible metal to insulator transition 
[loj . At very high pressures compression is sufficiently 
strong that core electrons start to overlap and Pauli ex- 
clusion principle favours less symmetric charge distribu- 
tions with p-d hybridizations [a, 0] . The Coulomb repul- 
sion between core and valence electrons leads to charge 
localization in interstitial positions, instead of around nu- 
clei, and finally to an insulating state From a sim- 
ple, uncorrelated, reflective metal, sodium turns into a 
complex, correlated, transparent insulator. It is a whole 
change of paradigm [ll|. 

The crystal structure of the recently discovered insu- 
lating state was identified as a distorted, six-coordinated 
double- hexagonal closed-packed structure (hP4) . Ex- 
perimentally, it is found that under pressure body- 
centred cubic Na first becomes face-centred cubic (at 65 
GPa) and then transforms into a cI16 structure [3] (at 
103 GPa). Besides adopting many different structures 



near the minimum of the melting curve [Sl (118 GPa), 
it crystallizes in the oP8 and tI19 phases [6( (above 125 
GPa), before undergoing a final transition to the newly 
proposed hP4 phase at around 200 GPa. Theoretically 
5], the incommensurate tllQ and the oP8 phases were 
found to have nearly equal enthalpies and the transition 
to the hP4 structure occurs at slightly higher pressures 
(260 GPa). Measurements were performed at room tem- 
perature, while the zero-temperature structural calcula- 
tions neglected entropic contributions and lattice dynam- 
ics (phonon enthalpies along with quantum effects on the 
light sodium ions), and hence could not provide a direct 
confirmation of the theoretical results [12]. 

In the present Letter, we study the optical response 
of the hP4 phase employing the highly accurate Bethe- 
Salpeter equation (BSE), which is the state-of-the-art ap- 
proach for first-principles predictions of optical proper- 
ties in solids |13j . Our results provide evidence for the 
emergence of relevant excitonic (i.e. electron-hole corre- 
lation) effects, with the formation of an unusual charge- 
transfer exciton, and a strong anisotropy in the dielectric 
response. They represent a stringent fingerprint of the 
hP4 crystal structure of transparent dense sodium. 

Absorption spectra and excitons - The BSE can be cast 
into an effective two-particle Schrodinger equation for the 
wavefunction x{Yh,Ve) of the electron- hole pair (i.e. the 
exciton) : 



Hp 



(1) 



where the excitonic Hamiltonian H^xc includes the 
electron- hole interaction [13]. In our first-principles ap- 
proach, is represented on a Kohn-Sham basis: 



vl';,(r,„r,) = ^vl,k-0:k(r^)</>ck(re), (2) 
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where v (c) runs on valence (conduction) bands and k 
is in the first Brillouin zone. The optical spectrum is 
then obtained from the imaginary part of the dielectric 
function €2 = Im e: 



e2{uj) = lim ^ V I 



(3) 

When excitonic effects are unimportant the eigenval- 
ues E\ correspond to vertical transition between one- 
quasiparticle energies ijJek ~ -S^ik (which we calculate in 



the GW approximation [ij]) and the coefficients 



become diagonal, reducing Eq. ([3|) to the well known 
Fermi's golden rule Q. We investigate the effect of pres- 
sure on the formation of excitons, through charge local- 
ization and the reduction of the screening of the electron- 
hole interaction. Excitons can be identified by means of 
the twofold role that they play on the absorption spectra. 
First, the modification of the excitation energies induces 
a shift of the positions of the peaks in the spectrum: as a 
consequence, new peaks can appear inside the fundamen- 
tal band gap (bound excitons) ; second, the mixing of the 
independent-particle transitions leads to a global reshap- 
ing of the spectrum towards lower energies (continuum 
excitons). Our ah initio solution of the BSE proceeds 
in three steps. We start from a ground-state calculation 
using density-functional theory in local-density approxi- 
mation (LDA), where we employ the lattice constants of 
the hP4 phase from Ref. la corresponding to representa- 
tive pressures (see Ref. [l5| Tab. 1). We treat 3s, 2p and 
2s electrons explicitly as valence, leaving only the very 
tightly bound Is in the frozen ionic core (this required 
us to use a energy cut-off of 80 Ha in the plane-wave 
representation of the wave functions) . We then calculate 
the quasiparticle GW electronic structure that we use as 
input for the Bethe-Salpeter equation ([T]). In the calcu- 
lation of the screening and the GW corrections, where 
we adopt the Godby-Needs plasmon-pole model [16|, we 
use 150 bands and a cut-off of 20 Ha on the wave func- 
tions. To converge the absorption spectra up to 7-8 eV 
we consider 2 valence and 8 conduction bands and sample 
the first Brillouin zone with a uniform grid of inequiva- 
lent 2000 k points, slightly shifted from high-symmetry 
k points. We then apply a 0.1 eV gaussian broadening 
24|. 

Excitons in dense sodium.- To illustrate our discus- 
sion, we first consider the results obtained with the lattice 
parameters calculated by Ma et al. [1] in correspondence 
to a pressure of 320 GPa. We then analyse the depen- 
dence of the excitonic effects with applied pressure. In 
Fig. [TJa) we compare the absorption spectra that we 
obtain by solving the BSE or by using Fermi's golden 
rule (averaged with respect to the different light polar- 
izations). In the latter case transition energies are calcu- 
lated either using LDA or the GW approximation. The 
GW corrections on the band structure induce an almost 




FIG. 1: (Color online) (a) Excitons in BSE (red curve, main 
panel) change substantially the shape of the GW (blue) and 
LDA (green) spectra calculated without electron-hole inter- 
actions. They redshift the absorption onset (grey arrow) by 
~0.1 eV with respect to the fundamental GW quasiparticle 
gap. The apparent shift of the spectrum is only partially 
accounted for by the change of the JDOS, from the GW ex- 
citation energies to the density of excitonic states of the BSE 
solution (inset). The remaining effect is due to mixing of 
formerly independent valence-conduction transitions, (b)-(c) 
Electron-hole pair wavefunction ^A(r?i,re) for the first exci- 
ton. When the hole (black dot) is located either near a corner 
atom (b) or in interstitial position (c), the electron charge is 
always on the central atoms (isocontour in red: plot of |^a|^ 
at fixed at 0.3 of its maximum), (d)-(e) Kohn-Sham one- 
particle real wave functions 0„k(r), top valence (d) and bot- 
tom conduction (e) states at the F point: isocontours at -1-0.3 
(red) and —0.3 (blue) of their maximum. 



rigid blue-shift to the LDA spectrum, reflecting the cor- 
rection of the usual underestimation of the quasiparticle 
band gap in LDA (see Ref. 15| Fig. 1). The onset moves 
from 1.7 eV in LDA to 3.0 eV in GW. It corresponds 
to a direct transition in F between the top valence and 
the bottom conduction bands. Adding interaction be- 
tween GW-quasiparticles, i.e. quasiparticle electron-hole 
attraction, we observe a substantial modification of the 
spectrum: i) strong redistribution of the optical oscilla- 
tor strength to lower energies, changing the global shape 
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of the spectrum; ii) strong narrowing of the main peak at 
~ 5 eV and iii) appearance of bound excitons inside the 
gap. It turns out that electron-hole attraction reduces 
the transition energies by about 0.1 eV, determining a 
global red-shift of the joint density of states (JDOS) (in- 
set of Fig. lUa)). The absorption onset is now at 2.9 eV 
(grey arrow in Fig. [Ija)). This implies the formation 
of a bound exciton in the fundamental gap with a bind- 
ing energy of about 100 meV, that is larger than the one 



13|. 



obtained in standard semiconductors 

In Fig. [IKb)-(c) we analyse the shape of the wave func- 
tion corresponding to this bound electron-hole pair. We 
fix the position of the hole in a particular position and 
determine the corresponding density distribution of the 
electron. Our calculation reveals that this bound exciton 
is, in particular, a charge-transfer exciton. While the hole 
is created in the vicinity of one corner atom (Fig. [Hb)), 
the excitation has largely moved the electron charge to lo- 
calize on the central atoms, forrning a separate sub-lattice 
over many unit cells (see Ref. [la] Fig. 2). Remarkably, 
the same electronic distribution is obtained even when 
the hole is located in the interstitial regions (Fig. HJc)), 
where the valence electrons accumulate when they are 
forced away from the atoms 0, [llj . Charge-transfer ex- 
citons are typical of ionic insulators like alkali-halide crys- 
tals [l3 . They are intermediate between strongly bound 
Frenkel and weakly bound Wannier excitons, as they in- 
volve the transfer of one electron from a negative ion to 
a positive ion. In the present case, in particular, the role 
of the anion is played by the interstitially confined elec- 
tron gas. In this sense, transparent sodium behaves as 
an unconventional inorganic electride [ISJ . characterised 
by the simplest possible anion: an electronic charge with- 
out ionic core. The strong enhancement of the peak at 
around 5 eV, which is the mark of the formation of a 
resonant exciton, instead is not explained by the simple 
modification of the transition energies, since the JDOS 
curves remain very similar. It is rather due to the mix- 
ing of a very large number of independent electron-hole 
transitions in the same energy range (see Ref. Fig. 
1), each of them participating with a small contribution 
to the excitonic peaks. It is the sign of highly correlated 
excitonic states. 

Anisotropy and inhomogeneity.- By comparing the op- 
tical spectra obtained with light polarized either along 
the hexagonal c axis or in the perpendicular ab plane 
(Fig. [21(a)), we find that the spectra with the two light 
polarizations remain very different, both with and with- 
out the inclusion of excitonic effects. In particular, in 
the in-plane direction, sodium is completely transparent 
up to 4.5 eV. In fact, the absorption onset and the first 
structures in the spectrum in Fig. [Ija) are entirely due 
to absorption of light polarized along the c axis. Also 
the exciton represented in Fig. [IKb)-lUc) is dark in the 
in-plane direction. This behaviour can be understood 
as due to the particular symmetry of the one-particle 
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FIG. 2: (Color online) (a) Spectra calculated with polariza- 
tion either along the hexagonal c axis (blue) or in the in-plane 
ab directions (red). In both cases, dotted lines are for GW 
spectra without electron-hole interaction and solid lines for 
BSE spectra with electron-hole interaction, (b) Absorption 
spectra calculated at lattice parameters corresponding to in- 
creasing pressures (data from Ref. [^), with or without exci- 
tonic effects (solid and dotted lines, respectively). 



wave functions involved in the transitions. For instance, 
by symmetry at the F point the transition between the 
top valence (Fig. [T](d)) and the bottom conduction (Fig. 
[IJe)) has zero oscillator strength in the in-plane direction. 
On the other hand, inhomogeneities in the charge distri- 
butions are responsible for the so-called crystal local-field 
effects [l^. In general, an external field can induce spa- 
tial charge fluctuations that are rapidly varying on the 
microscopic scale. It is intuitively clear that the self- 
consistent response to these induced local fields becomes 
more important the more the system is inhomogeneous 
and polarizable. Whereas in the absorption spectra of 
bulk solids local fields are generally weak [l^], in dense 
sodium they induce a sizeable effect (see Ref. [l^] Fig. 
3), particularly for the main peak and c axis polariza- 
tion (while the onset is unaffected). This is a further 
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indication that the charge distribution has become very 
inhomogeneous . 

Dependence on the pressure - The apphcation of in- 
creasing pressure induces as a first effect a band gap 
opening in the quasiparticle band structure [5j]. As a 
consequence also the optical spectra move to higher en- 
ergies (Fig. mb)). Moreover, the quasiparticle gap from 



indirect becomes direct at the F point (see Ref. [15 1 
Fig. 5). Additional compression leads to localize more 
the electronic charge and to reduce the screening of the 
interaction between the electron and the hole. Both con- 
tribute to increase the excitonic effects: the electron and 
the hole on average stay closer and their effective interac- 
tion is stronger. This explains why the largest differences 
between the spectra with and without electron-hole in- 
teraction occur for the spectrum at the highest pressure 
(Fig. Eljb)). In fact, the static dielectric constant, which 
we calculate in random-phase approximation (RPA) [l3| , 
decreases with increasing pressure, passing from 18 at 
220 GPa to 7 at 420 GPa. In addition, we consider the 
dynamical screening, — Ime~^(q, w), which can be mea- 
sured by electron-energy loss spectroscopy (EELS) or in- 
elastic x ray scattering (IXS) and which we calculated 
in time-dependent LDA (2^ (with the inclusion of 40 
bands) . The EELS spectrum, for each momentum trans- 
fer q (either along the c axis or in the perpendicular ab 
plane), is characterised by a main structure, which cor- 
responds to the frequency of the collective oscillations of 
the valence electrons in the system, i.e. plasmons (Fig. 
[3]). These plasmon structures show a small dispersion 
with the momentum transfer q (see Ref. [ll] Fig. 6). 
While detailed features depend on band structure prop- 
erties, in this case the Drude model is still partially able 
to track the dependence of the plasmon frequency ujp 
with the compression. According to the Drude model, 
the plasmon frequency increases linearly with the square 
root of the density [l| : Wp = ^/A^^p. Using this formula at 
the densities corresponding to the three considered pres- 
sures we have Wp = 13.5, 14.6 and 15.4 eV (grey arrows 
in Fig. [3]) when we consider only 4 valence electrons in 
the unit cell. In a way, this corresponds to a definition of 
"valence" electrons in this extreme situation where the 
usual view of electrons belonging to a particular atom is 
no more strictly valid. Being based on a uniform electron 
gas, the Drude model describes an average behaviour of 
the dielectric response of the system, but is obviously not 
able to catch its anisotropy, which leads to the splitting 
between the plasmons along the c axis and in the perpen- 
dicular plane (the Drude plasmon energy falls somewhat 
in the middle). Moreover, at low pressures a smaller 
further plasmonic structure appears only in the in-plane 
direction at about 8 eV, which is strongly suppressed 
by increasing the applied pressure (see Ref. [l5| Fig. 
7). Both the anisotropy in the lattice structure and the 
higher-energy localisation of the spectral weight for the 
in-plane polarization (see Ref. flB] Fig. 7) contribute in 




FIG. 3: (Color online) EELS spectra, — Im (q, u), along 
the two polarizations for a momentum transfer q equal to 
0.1 reciprocal lattice units at 3 different pressures: 220, 320, 
and 420 GPa. The grey arrows mark the plasmon energies 
obtained using a Drude model for a homogeneous electron 
gas with a density corresponding to 4 valence electrons per 
unit cell. 



the same way to the splitting of the plasmon energies in 
the two directions. 

In conclusion, our first-principle predictions for both 
absorption and EELS spectra can be indeed used to con- 
firm the crystal structure of this phase of dense transpar- 
ent sodium. First of all, since the binding energy of the 
bound exciton is of the order of a hundred of meV, our re- 
sults for the optical spectra support the hypothesis of Ma 
et al. , which was based only on ground-state and band 
gap calculations. Moreover, even though ground-state 
calculations overestimate the pressure value at which the 
transition to the hP4 phase occurs our predictions 
for the optical properties are not biased by this error. In 
fact, while excitonic effects are varying with the applied 
pressure (and can be influenced by the temperature as 
well [llj), from our analysis it turns out that key fea- 
tures of the spectra like the polarization anisotropy are 
essentially linked to the symmetry of the crystal struc- 
ture and can be observed in a very large range of pressure 
(see Ref. 



15| Fig. 4). This also implies that sodium be- 



comes transparent in the visible at a lower pressure in 
one polarization direction while being still "metallic" in 
the other. Therefore we suggest the measurement of the 
spectra along the two polarizations as the key test to 
verify conclusively that transparent sodium crystallizes 
in the hP4 structure. Photoluminescence would put in 
evidence the nature of the bound exciton predicted in 
the present work and that has a low oscillator strength 
in optics. 
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P 


2.874 


4.238 


220 


2.838 


4.057 


260 


2.769 


3.817 


340 


2.705 


3.668 


420 



TABLE I: The lattice parameters of the Na-hP4 phase used in the present work (in A) at different pressures (in GPa), from 
Y. Ma et al., Nature 458, 182 (2009). 




FIG. 1: Band structure at 320 GPa. The GW corrections (blue) induce an almost rigid upshift of the conduction states 
calculated in LDA (red). The zero of the energy scale is set at the top- valence energy. The fundamental gap is direct at F: 1.7 
eV in LDA and 3.0 eV in GW. In yellow: region of the Brillouin zone where the vertical transitions from top valence to bottom 
conduction states have zero dipole oscillator strength in the in-plane direction. They contribute only to the absorption for light 
polarized along the hexagonal c axis, forming the first peak in the spectrum at lower energies. In light-blue: zone where they 
contribute to both polarizations (in this case the top valence band is twofold degenerate). 
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FIG. 2: Exciton wavefunction. Isocontour plot of the exciton wavefunction re) at 0.3 of its maximum. The position 

Th of the hole is fixed on one corner atom (black dot). The isocontour represents the electron charge extending over many unit 
cells (here shown a 3x1x3 supercell). This implies that the electron and the hole are not strongly bound (as it would be for a 
Frenkel exciton). Na behaves as a "bipartite" system: the hole lies on a one site and the electrons on the others. 
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FIG. 3: Crystal local-field effects. Dotted lines: absorption spectra calculated with Fermi's golden rule using GW quasipar- 
ticle energies. Dot-dashed lines: with the inclusion of crystal local fields. The effect is quite strong, in particular for the c axis 
polarization (blue). The spectral weight is blueshifted (but the absorption onset does not change) and its intensity reduced. 
In sp semiconductors the effect of local fields on absorption spectra normally is less relevant than what found here. Local 
fields become more important when the system is less homogeneous (but still polarizable) . This is a manifestation of the very 
inhomogeneous distribution of the electronic charge. 
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FIG. 4: Polarization dependence of the transparency. Yellow zone: energy range of visible light (up to 3.1 eV). At low 
pressures Na is transparent in the in-plane direction (solid lines) and reflective along c axis (dot-dashed lines). By increasing 
pressure, the peak in the in-plane direction (solid lines) is quite rigidly blueshifted and its intensity is reduced accordingly (as 
an effect of the /-sum rule). Instead, for the polarization along the c axis (dot-dashed lines) at low pressures the spectrum 
is characterised by a double peak, while by increasing pressure the first peak progressively looses strength developing as a 
shoulder of the second one. 
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FIG. 5: From indirect to direct semiconductor. Evolution of the LDA band structure with appHed pressure. At lower 
pressures Na has an indirect gap (the bottom conduction is at T and the top valence along T — K and F — M). By increasing 
pressure, the gap becomes direct at F. The change occurs at the pressure P=260 GPa where in LDA the hP4 structure becomes 
energetically more stable than the tI19 one. While the valence bandwidth is quite insensitive to pressure, there appears a clear 
band gap opening. The zero of the energy scale is located in correspondence to the top valence energy. 
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(o[eV] 

FIG. 6: Plasmon dispersion. EELS spectra calculated in TDLDA for increasing momentum transfers q, eitiier along the c 
axis (blue) or in the perpendicular ab plane (red). The lattice parameters correspond to a pressure of 320 GPa. In the direction 
of the arrow, at each curve the momentum transfer is increased by 0.1 reciprocal lattice units, from 0.1 to 0.5. Both the plasmon 
energy and the intensity of the plasmon peaks reduce with increasing q, but the dispersion with q is not very large. The double 
plasmon peak in the ab plane (red) tends to broaden into a single structure. This means that the polarization anisotropy of 
the loss function is robust with respect to the momentum transfer. 
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FIG. 7: Dielectric function. Real and imaginary parts of the dielectric function e = ei + ie2 for in-plane polarization (a) at 
220 GPa and (b) 420 GPa, and c axis polarization at 420 GPa (c) for a momentum transfer q = 0.1 reciprocal lattice units. 
The grey arrows mark the plasmon peaks in the energy loss spectra L = e2/{e\+e\). They correspond to zeros of e\ where also 
£2 is small. The first plasmon at 7 eV at the lowest pressure (a) is strongly suppressed by increasing the applied pressure: £2 
is rising in intensity in the same energy region (b), as an effect of band gap opening. The plasmon peak in the c axis direction 
is located at lower energies. In fact, in the in-plane direction the main £2 peak has a shoulder extending to higher energies (b) 
that is absent in the c axis direction (c). 



